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ABSTRACT 

We report new observations of the fundamental J = 1 — transition of HC1 
(at 625.918GHz) toward a sample of 25 galactic star-forming regions, molecular 
clouds, and evolved stars, carried out using the Caltech Submillimeter Observa- 
tory. Fourteen sources in the sample are also observed in the corresponding H 37 C1 
J = 1 — transition (at 624.978GHz). We have obtained clear detections in all 
but four of the targets, often in emission. Absorptions against bright background 
continuum sources are also seen in nine cases, usually involving a delicate balance 
between emission and absorption features. From RADEX modeling, we derive 
gas densities and HC1 column densities for sources with HC1 emission. HC1 is 
found in a wide range of environments, with gas densities ranging from 10 5 to 
10 7 cm -3 . The HC1 abundance relative to H 2 is in the range of (3 - 30) x 1(T 10 . 
Comparing with the chlorine abundance in the solar neighborhood, this corre- 
sponds to a chlorine depletion factor of up to ~400, assuming that HC1 accounts 
for one third of the total chlorine in the gas phase. The [ 35 C1]/[ 37 C1] isotopic ratio 
is rather varied, from unity to ~5, mostly lower than the terrestrial value of 3.1. 
Such variation is highly localized, and could be generated by the nucleosynthesis 
in supernovae, which predicts a 37 C1 deficiency in most models. The lower ratios 
seen in W3IRS4 and W3IRS5 likely confine the progenitors of the supernovae to 
stars with relatively large mass (>25M Q ) and high metallicity (Z~0.02). 

Subject headings: ISM: abundance - ISM: molecules - submillimeter: ISM 
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Introduction 



Initially of interest as a potential coolant for molecular clouds, hydrogen chloride (HQ) 



was first detected in the J = 1—0 transition at 625.918 GHz 



De Lucia. Helminger. Gordv 



Blak e. Keene. fc Phillips ( 119851 ). and in 



Zmuidzinas et al. 



( 119951 ). both aboard the 



19711 ) toward IRc2 in Orion A, in emission, by 
absorption against dust continuum in Sgr B2, by 
NASA Kuiper Airborne Observatory with a beam width of 2(3. Ground-based observations 
at a higher spatial resolution of 12" were carried out using the 10.4-m telescope of 



Schilke. Phi 



lips, fc Wand (119951) and 



Salez. Frerking. fc Langer 



the Caltech Submillimeter Obse rvatory (CSO) by 
Salez Frerking. fc Langerl fll996f ) toward OMC-1 and MonR2. 
(119961 ) also detected the H 37 C1 J = 1 — transition toward OMC-1. These studies found the 
HC1 abundance to be 10~ 9 — 10~ 10 relative to molecular hydrogen, indicating that chlorine 
is depleted by a factor of several hundred compared to its sola r neighborhood abundance , 



Salez. Frerking. fc Langer 



and t hat HC1 is not a significant coolant of molecular clouds. 
( 119961 ) also reported a [ 35 C1]/[ 37 C1] isotopic ratio that is 1.3 — 2 times higher than the solar 
neighborhood value, indicating variations of the isotropic ratio in the Galactic interstellar 
medium. 



Chemistry of chlorine and 
tigat e d by a number of authors 



lydrqgen chloride in the interste 



1986 



Schilke. Phillips, fc Wang 



Jura 



1995 



1974 



Dalgarno et al. 



Neufeld fc Wolfire 



1974 



lar medium has been inves- 



Blake. Anicich. fc Huntress 



20091 ). establishing a relatively 



simple chemical network involving less than a dozen other atomic and molecular species, 
as well as electrons and cosmic-ray particles, in the production and destruction of HC1. In 



dense molecular clouds, such a network yields, among o ther species, predominant 



taking up about one thi rd of the chlorine in gas phase ( jSchilke. Phillips, fc Wang 



Neufeld k Wolfire 



200% . 



y HC 1 



1995 



HC1 J 



transition falls on the shoulder of a relatively strong atmospheric water 
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absorption line. Consequently, ground-based observations of HC1 have been limited to the 
rare occasions of superb weather at high sites, such as Mauna Kea, and relatively little is 
known about HC1 in different categories of objects, such as evolved stars, active star forming 
regions, molecular clouds with no apparent star formation activity, etc. In particular, the 
[ 35 C1]/[ 37 C1] isotopic ratio is of interest in different categories of objects, as variations of 
this ratio can be produced through nucleosynthesis in massive stars followed by mixing o f 



the p rocessed material into the surrounding molecular cloud (jSalez. Frerking. fc Langer 



19961 ). The chlorine isotopic ratio could thus offer insights into the evolutionary history of 



the molecular material. 

To address these questions, we took advantage of several periods of superb weather on 
Mauna Kea to conduct a comprehensive survey of HC1 in a sample of well known objects. 
In the following sections, we first describe the observational details, then present the results 
and data analysis, followed by discussions and conclusions. 



2. Observations 

Observations described here were carried out on September 11, 2008, March 2, and 
April 7-15, 2010 with the 10.4 meter CSO telescope. The 690 GHz facility heterodyne 
receiver was tuned to 625.918756 GHz and 624.977821 GHz, respectively, for observations 
of HC1 J = 1 — and H 37 C1 J = 1 — lines. The lines fall on the shoulder of a prominent 
atmospheric water absorption line, therefore requiring premium weather conditions. 
The atmospheric opacity measured at 225 GHz was mostly 0.03-0.04, and the typical 
double-side-band system temperature was 2500 K. There were marginal periods in April 
7-15, 2010 when atmospheric opacity was around 0.05-0.06, resulting in a typical system 
temperature of 6000 K. Both the FFTS and 4 GHz AOS facility spectrometers were used in 
1 GHz mode as backends. 
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The Dish Surface Optimization System (DSOS) (ILeong et al.ll2006l ) was used to correct 
gravitational deformations of the telescope surface, maintaining a clean 13'.'5 beam over a 
wide range of elevations. Temperature calibration using a hot chopper wheel yields TjJ, 
the antenna temperature corrected for atmosphere and telescope hot spillover. A main 
beam efficiency of ~50% measured on Jupiter and Uranus is used to convert T* A to T m b, the 
main beam brightness temperature. Telescope pointing was checked on Jupiter, Mars, and 
Uranus, and was found to be within 4". 

Observations were conducted in the beam switching mode at a frequency of 0.523 Hz 
and with a chop throw of 300". Typical integration time on each source was 300 seconds, 
resulting in a rms of ~0.15 K at a velocity resolution of 0.31 km s -1 . The source selection 
was primarily based on its accessibility at reasonable elevations, with a secondary goal 
of covering a diverse collection of objects. Due to its high critical density (~10 8 cm -3 ), 
the excitation of the HC1 J = 1 — transition is most likely sub-thermal, and radiative 
excitation, in the presence of continuum emission, is likely to play an important role in the 
excitation of the transition. The presence of dense gas and/or continuum radiation in an 
object would likely yield a detection. In the end, 25 sources/positions were observed, among 
them active star forming regions, giant molecular clouds with no apparent star forming 
activities, evolved stars, shocked gas in a SNR, and ultra-compact HII regions. Table [H lists 
observational details for all observed sources. 



3. Results 

For sources with matching HC1 and H 37 C1 J = 1 — spectra, we present in Figures [IH3] 
HC1 spectra in the left column and H 37 C1 spectra in the right column. The antenna 
temperature scale is kept the same for each source for easy comparison of the line intensities. 
The velocity scale, Ulsr ~~ Vo, refers to the main hyperfine component; values of Vq are 
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denoted in each panel. Toward IRc2 in OMC-1, the HC1 and H 37 C1 line profiles are blended 
heavily with spectral lines from other species, in the case for H 37 C1 predominantly methanol. 
It is not easy to untangle the blending and derive useful information from the resulting line 
profile. Instead we opted to observe two positions south of IRc2 with declination offsets 
of -60" and -100", respectively. The two pos i tions correspond to local maxima in the HC1 



with the CS clump LSI ( IMundy et al. 



map of OMC-1 by ISchilke. Phillips, fc Wangj (119951 ) . The (0, -60") position also coincides 



19861 ). and the (0, -100") position Orion S. Figure 1 



show the HC1 J = 1 — spectra for the remaining sources in the sample. 



3.1. Hyperfine Structure Fitting 

The HC1 hyperfine structure, resulting from the interaction between the electric field 
and chlorine nuclear spin, splits the J = 1 — transition into three components, with 
the outer two separated by —6.35 and 8.22 km s -1 respectively from the strongest middle 
component for H 35 C1. For H 37 C1, the separations are —5.05 km s -1 and 6.45 km s _1 , 
respectively. Statistical weights for the upper level of the hyperfine transitions dictate 
optically thin line strength ratio of 1:3:2. In the observed sample, the three-component 
hyperfine pattern is clearly seen in many sources, but the components often do not conform 
to the optically thin ratio, indicating deviations from the optically thin limit or differences 
in excitation between the three components. 

Under the assumption that the hyperfine components share the same excitation and 
kinematic characteristics, fits to t he line profiles can ind ependently provide the optical 



depth and excitation temperature 



Kwan &: Scoville 



19751 ). Such a utility is readily available 



in CLASS, the spectral line data reduction package developed and maintained by IRAM. 
The HFS fit in CLASS provides four parameters: AT ■ r main , Vlsr, w the line width, 
and r main . r main is the total optical depth from all hyperfine components of the line, and 
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AT = T ex — Tb g is the difference between the line excitation temperature and background 
temperature. Table |2] lists results of such hyperfine structure fits for sources in our sample 
that show the J = 1 — HC1 or H 37 C1 emission. Cases involving absorption are more 
complex, as seen in G34.26+0.15 and W43-MM1 in Fig-|3]for example. They are often 
the result of a delicate balance between emission and absorption features. We defer the 



treat ment and discussions for sources with absorption to a separate paper ( lYoshida et al. 



201fl). 



In fitting the hyperfine profiles in T m b scale, one needs to be mindful that the line 
emission may suffer from beam dilution when a source is smaller than the telescope beam. 
The uncorrected line temperatures would affect the excitation temperature and optical 
depth derived from HFS fit of the HC1 line profile. Thermal dust emission as seen in 
submillimeter continuum maps provide a good measure for gas distribution in the source. 
We compare the size of the 350^tm dust cores to the telescope beam of 13'.'5 to determine 
whether the HC1 emission from a given source suffers appreciable beam dilution. Most of 

I — ^^Pn 

our sources in the sample are observed by iDotson et al.l (120101 ) with the SHARP camera at 
the CSO, and the FWHM of the dust cores range from ~30" to ~60". We therefore do not 
consider beam dilution a significant source of error for our sample of objects in deriving the 
excitation temperature and optical depth for the HC1 J = 1 — transition. One exception 
may be G5.89-0.39, where SMA continuum image at 870/zm reveals a dust core ~ 6" — 8" 



in size (ITang et al. 



20091 ). implying a factor of ~ 3 in beam dilution in the observed line 



temperature. 

The HFS fits are sensitive to the signal-to-noise ratio in the observed line profiles. 
AT and hence T ma i n are not well constrained for relatively noisy line profiles. For sources 
with both HC1 and H 37 C1 detections, one would expect that the values of AT derived from 
HFS fits of the two transitions be very close. In reality, they tend to agree within error 
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bars, albeit rather large error bars in some cases, resulting in questionable estimates of 
the optical depths. We attempted to derive alternative estimates of the optical depths for 
both HC1 and H 37 C1 lines by imposing a fixed value of (AT)*, a weighted average of the 
values obtained for the two lines. Using (AT)* and the T m b from Gaussian fits of the three 
components in each profile, we derive r(HCl)* and r(H 37 Cl)* as listed in Table EJ The error 
of T m b/(AT) amplified by e T translates into the error of r* for each components, resulting 
in relatively large errors for r* derived from this approach when individual components are 
optically thick. We nonetheless consider r* a more reliable measure of the optical depth 
than that obtained directly from HFS fits. For sources that have only HC1 observations, 
(AT)* and r(HCl)* in Table [2] are from direct HFS fits. For G5.89-0.39, we also present 
fits corrected for the likely case of beam dilution. Compared to the uncorrected case, the 
differences are mostly in the line excitation temperatures, with optical depth remain largely 
unchanged. 

IRC+10216 shows a rather distinctive line profile. It is likely the result of convolution 
of the hyperfine structure with a shell profile often seen in circumstellar envelopes. As an 
approximation, we chose to fit two Gaussian profiles for IRC+10216 at LSR velocities of 
—17.8 km s _1 and —27.4 km s _1 , respectively, to simulate the shell profile in the hyperfine 
fits. As shown by the red curve in the appropriate panel in Figure El this approach 
adequately fits the line profile, yielding a solution that indicates that the HC1 emission 
originated from a circumstellar shell with an expansion velocit y of 4.8 km s -1 , compared to 
an observed terminal velocity of 14.5 km s _1 (see for example 



Schoier fc Olofsson 



d20QlO). 



The [H 35 C1]/[H 37 C1] isotopic ratio in the shell is likely in excess of 4 from the optical depth 
shown in Table |2j 

CIT6 is widely regarded as similar to IRC+10216 with com parable chemical and 



physical characteristics, but a lower dust-to-molecular mass ratio ( IZhang. Kwok. fe Trung 
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20091 ). HC1 is no exception in that regard. It is clearly detected but the line is too weak 
to allow for a good HFS fit. Nevertheless, this observation serves as an additional proof 
that HC1 as a species exists in detectable quantities in circumstellar shells of evolved 
stars. IC443G1 is a dense core in a ring of shocked material associated with SNR IC443 



( Ivan Dishoeck. Jansen. fc Phillips 



19931 ). There are tantalizing signs of HC1 detection, but 



the line is again too weak to derive any useful parameters. 

G5. 89-0. 39 is another interesting case with a rather distinctive line profile. Both 
HC1 and H 37 C1 line show a pronounced red wing, indicating presence of higher velocity 
components. A three-component fit, as shown by the red curve in the top left panel of 
Figure |3j represents an excellent match to the observed line profile. As listed in Table |2j 
the three components are at LSR velocities of 9.3, 17.9, and 38.3 km s" 1 , with FWHM 
line widths o f 4.9, 13.5, and 5.4 km s _1 , respectively. From the interferometric study of 



Hunter et al. 



( 120081 ). the broad (~13 km s _1 ) line emission is seen in C 34 S, SiO, CH 3 OH, 
HC3N, C 17 0, and H 13 CO + toward the dust emission source SMA-N, in the north-northeast 
extension of the molecular ridge that bounds the northern and western edge of the UCHII. 
The 9.3 km s _1 feature likely originates from the northwest and southwest section of the 
molecular ridge. While the origin of the 38.3 km s _1 feature is not immediately clear, 
velocity features as highly redshifted as 78 km s _1 and blue-shifted as -61 km s _1 , mostly 
from various masers, were quoted by these authors in establishing a north-south bipolar 
outflow. In deriving the alternate optical depth of these velocity components in Table El we 
opted to use the AT from the HC1 HFS fits for the 17.9 km s _1 and 38.3 km s _1 components 
due to their low S/N ratio in the H 37 C1 line profile. From the derived r* in Table [2], the 
[HC1/H 37 C1] isotopic ratio in the source is most likely around 2, though the broad emission 
feature may have a higher ratio on the order of 3-4. 

Sgr A*, the dynamical center of the Galaxy, yields no HC1 detection to an rms level of 
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0.1 K at a resolution of 0.31 km s -1 . Similarly, Gl. 6-0. 025, a giant molecular cloud in the 
Galactic center region with no apparent sign of star formation, also yields a non-detection. 
Sgr B2(M) and Sgr B2(N), on the other hand, show broad strong absorption at both HC1 
and H 35 C1 J = 1 — transitions. 

From Table [21 it it clear that the [H 35 C1]/[H 37 C1] isotopic ratio, based on that of the 
optical depth, is rather varied, from 4 and over in Mon R2, IRC+10216, and DR21(OH), 
to around 3 in OMC-1, W3(OH), and IRAS16293, down to around 2 in G5. 89-0.39, and to 
nearly reaching parity in W3IRS4 and W3IRS5. In comparison, the terrestrial [ 35 C1]/[ 37 C1] 
abundance ratio is ~3.1. 



3.2. RADEX Simulations 



In addition to the optical depth obtained from hyperfme structure fitting, we can derive 
other physical parameters from the line profiles by running LVG simulations to fit the three 
HC1 hyperfme componen ts. RADEX, a computer program developed and maintained by 



van der Tak et al. 



(120071 ) for fast non-LTE analysis of interstellar line spectra, was used in 
the LVG mode to find the best combination of gas density and HC1 column density that 
could reproduce the observed line temperatures under given gas temperature. 

The molecular data file for HC1 under RADEX includes the hyperfme levels of the 
lowest 8 rotational levels, up to 839 K above ground level, with rates for radiative and 
collisional transitions among t hem. The collisional rat es between HC1 and H 2 are scaled 



from those of HC1 and helium ( Neufeld fc Green 



1994T ) . The average Galactic background 



radiation field adopted in RADEX consists of the cosmic microwave background, average 
Galactic starlight in the solar neighborhood in ultraviolet /visible/near- infrared part of 
the spectrum, and single-temperature fit to the Galactic thermal dust emission for the 
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far-infrared and submillimeter part of the spectrum. Contributions from non-thermal 
radiation in the Galaxy are also included for low frequency part of the spectrum (<30 
GHz). This composite interstellar radiation field and any additional continuum background 
from user input provide radiative excitation in the RADEX formulism. For the subset of 
our sample where HC1 emission is observed, no significant continuum is seen in any of 
these sources. We therefore use the default composite interstellar radiation field in RADEX 
simulations. 

Simulations were run in off-line mode to cover a large parameter space, and relative x 2 
was sought between the observed line temperatures and the model results relative to the 
observed line temperature. Under the assumption that the observed HC1 emission suffers 
no significant beam dilution, this process yields for a given gas temperature the most likely 
combination of gas density and HC1 column density to reproduce the observed line strength 
of the hyperfine components. Figure [5] shows a RADEX simulation run for OMC-1 (0, -60) 
with a gas temperature of 45 K. Contours are the relative \ 2 starting from 0.5 with a step 
of -0.1; blue and red curves mark the hyperfine component ratios under optically thin limit 
(3 between F=5/2-3/2 and F=l/2-3/2 transitions in blue, 1.5 between F=5/2-3/2 and 
F=3/2-3/2 transitions in red). The solid smooth line denotes the observed optical depth 
for the HC1 line and the two dotted lines its ±lcr error bar. In cases where the quoted 
errors for HC1 optical depth in Table [2] exceed 30%, we generally use 30% to set the error 
bar for the optical depth. The cross-section between the x 2 contours and the optical depth 
curves clearly confines gas density and HC1 column density in an often tight range. 

Table [3] lists for each source the adopted gas temperature and H 2 column density, 
the gas density and HC1 column density derived from RADEX simulations, as well as the 
HC1 abundance. Since there is no molecular data file for H 37 C1 to use under RADEX, we 
treat it the same as HC1 and fit the H 37 C1 line temperatures with the same procedure. 
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Instead of confining the ranges of gas density and H 37 C1 column density with the relative 
X 2 contours and optical depth curves, we use gas density ranges found in the corresponding 
HC1 simulation of the source in place of optical depth curves, ensuring that both HC1 and 
H 37 C1 emission originated from the same region. The resulting H 37 C1 column density is 
listed in Table [3] for sources with observed H 37 C1 line emission. For these sources we also 
list the 35 C1/ 37 C1 abundance ratio X[ 35 C1/ 37 C11 as measured by the ratio of HC1 and H 37 C1 



column densities. 



Cernicharo et al. 



(l2010bl ) recently suggested that the integrated intensity 
ratio be a good gauge for HC1 and H 37 C1 isotopic ratio. We listed in Table |3] the integrated 
intensity ratio, R, which appears to be in good agreement with our derived chlorine isotopic 
ratio. 



4. Discussion 

4.1. Where to Look for HC1 Emission? 

From Table [3] it is clear that HC1 emission is associated with warm dense gas with a 
density on the order of 10 6 cm -3 . Given a critical density of a few xlO 8 cm -3 , molecular 
excitation is mostly sub-thermal, dominated by radiative processes rather than collisions 
with H2. Sources with strong continuum and/or dust emission would likely yield strong 
HC1 line emission. 

The derived gas densities for the two OMC-1 cores appear to be somewhat higher than 



5 x 10 6 established by the multi-transitional study of HC3N (jBergin. Snell. &: Goldsmith 



19961 ). This may indicate that the gas temperatures adopted for these cores for the 
RADEX simulations should have been higher. Another factor may be the average Galactic 
background radiation field used in RADEX for radiative excitation. As HC1 transitions 
are short-ward of ~480/iin in wavelength, radiation from warm dust would contribute 
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significantly to the molecular excitation. Including local warm dust in the radiation field for 
RAD EX simulation would produce a more precise picture of the gas where HC1 emission 
originates. 



The non-detection toward Sgr A* li kely reflects a genuine 



circumstellar disk surrounding Sgr A*. 



Latvakoski et al. 



ack of HC1 molecules in the 



(119991 ) show a warm dust ring 



80" in size, as well as the "minispiral" within the 



ring. The dust 



ring coincides with a 



1983). With our 13'! 5 



molecular ring with sharp inner edges seen in HCN (IGiisten et al. 
telescope beam, we are clearly probing the molecular cavity around Sgr A*. Using a FWHM 
(15 km s _1 ) and the integrated line in tensity of the C 18 Q J= 2 — 1 line we obtain toward 



Sgr A*, and assuming LTE at 250 K ([Bradford et al. 



20051 ). we derive a C 18 column 



density of 9 x 10 15 cm 2 , and an H? column density of 5 x 10 22 with a C 18 abundance of 



1.7 x 10~ 7 ([Goldsmith. Bergin. fc Lis 



1997J) . A upper limit of 4 x 10 12 cm -2 is derived for 
the HC1 column density by RADEX simulations for an HC1 line of 0.1 K with a FWHM 
of 15 km s _1 (from the C 18 spectrum), yielding an upper limit of 8 x 10~ n for HC1 
abundance relative to H_2 toward Sgr A*, markedly lower than those listed in Table [3j 

The non-detection in G 1.6-0. 025 is intriguing. The source is a giant molecular cloud 
located in the Central Molecular Zone (CMZ) of the Galaxy, a region characterized 
by elevated gas temperature (typically ~70 K), higher gas den sity (> 10 4 cm" 3 ), and 



highly supersonic internal velocity dispersion (~15-50 k m s 1 ) (Morris & Serabyn 



Dense cloud cores are seen in various mo 



Gardner Boes 



1987 



Whiteoak k, Peng 



ecular tracers ( 



1989 



Gardner fc Whiteoak 



Kuiper et al 



1993 



199J). 



1981 



Menten et al. 



2009|), 



yet there is little sign of apparent star forming activity in the cloud. For the extended 
50 km s _1 cloud of Gl. 6-0. 025, we derive a C 18 column density of 2 x 10 15 from the 
C 18 Q J=2 — 1 spectr um taken at the CSO, assuming LTE at a gas temperature of 30 K 



( IMenten et al. 



20091 ). Again with a C 18 abundance of 1.7 x 10 7 , the cloud has an H 2 



-14- 



column density of 1 x 10 22 cm -2 . Similar to Sgr A*, we obtain an upper limit of 7 x 10 12 
for HC1 column density by modeling a HC1 line of 0.1 K with 28 km s _1 in line width, 
yielding an HC1 abundance of < 7 x 10~ 10 relative to H 2 . It would look as though the HC1 
non-detection in the 50 km s" 1 cloud of Gl. 6-0. 025 is the result of a combination of relative 
low H2 column density and HC1 abundance. 

It would be interesting to see if other clouds in CMZ have similarly low HC1 abundance 
seen in Sgr A* and the 50 km s _1 cloud of Gl. 6-0. 025. A prime target would be M-0.02-0.07, 
the so-called 20 km s™ 1 cloud ~2' northeast of Sgr A*. A comm on trait among the three 



clouds is the likely interaction with a SNR (IMenten et al 



20091 ). It is not clear though 



how such interaction would lead to relatively low production of HC1. 

In the same vein, IC443G1 may represent another case where chemistry plays an 



important role in the detection of HC 



material associated with SNR IC443 (Ivan Dishoeck. Jansen. fc Phillips 



Clump Gl is a dense core in a 



ring o f shocked 



1993J). It has 



a relatively rich chemistry compared to the pre-shock gas, but lacks many of the more 
complicated molecular species, such as methanol. These authors argued that such deficiency 
could be the result of disruption of chemical network brought on by the infusion of atomic 
hydrogen to the shocked gas by repeated shocks, as atomic hydrogen rapidly destroys many 
of the intermediate molecules in the reaction network. To derive an estimate of HC1 column 
density from the low signal-to-noise detection, we use RADEX simulation to reproduce an 



19931). 



HC1 line of 0.2 K with a FWHM of 10 km s" 1 (Ivan Dishoeck. Jansen. & Phillips 
The resulting HC1 column density of 3 x 10 12 cm" 2 , when compared to an H 2 column 



1993J), leads to an HC1 abundance of 



density of 1 x 10 22 (Ivan Dishoeck. Jansen. fc Phillips 
3 x 10~ 10 . Again, the weak HC1 line emission from IC443 Clump Gl may be the result of 
relatively low HC1 abundance and H 2 column density. More observations of similar cores are 
needed to confirm that shocked gas may have a somewhat lower HC1 abundance compared 
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to other similarly dense gas in star forming regions, PDRs, and Ultra-compact HII regions. 



Ultra-compact HII regions with reasonably high emission measures have proven to be 
a treasure trove for detecting HC1 absorption (see Figure Hj). Such cases typically involve 
a delic ate balance between a emission feature and an absorption feature. Our forthcoming 



paper (jYoshida et al. 



20101 ) will discuss these cases. 



4.2. The Case of IRC+10216 



In Figure |2] we shown that the HC1 line profile for IRC+10216 could be adequately fit 
by two Gaussians centered at —17.8 km s _1 and —27.4 km s -1 respectively, indicating that 
the HC1 emission is associated with a circumstellar shell with an expansion velocity of 4.8 
km s -1 . The HC1 emis sion is likely from a density-e nhanced shell at ~15 /; from the centra l 



star as seen in H CN (pchoier et al. 



and modeled by 



Corniner fc Millar 



2006|) 



2007) and SiO ( jSchoier. Olofsson. fc Lundgren 
( 120091 ). Radiation from the central star, as well as 
from local dust, contribute significantly to the molecular excitation. We therefore consider 
the HC1 abundance and [ 35 C1/ 37 C1] abundance ratio for IRC +10216 as listed in T a ble [3 



hig hly uncertain. Mo r e soph isticated modeling efforts such as 



and 



Wvrowski et a 



Cernicharo et al. 



Schoier fc Olofsson 



( 1200 lh 



(2006) are required to correctly interpret the results for IRC+10216. 
(l2010bf ) recently reported an HC1 abundance of 5 x 10~ 8 from their 



Herschel SPIRE and PACS observations of IRC+10216. It may be interesting to note that 
we derive an upper limit of 8 x 10~ 8 for the HCL abundance in CIT6 by treating the 
low signal-to-noise detection as a 0.2 K line with 20 km s _1 in line width. Similar to the 
procedure used in Sgr A* and Gl. 6-0. 025, the resulting HC1 column density is 5 x 10 12 cm -2 , 
with an adopted gas temperat ure of 40 K and an H 2 column density of 6.3 x 10 19 cm -2 



(jZhang. Kwok. fe Trung. 



200 
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4.3. 35 C1/ 37 C1 Abundance Ratio 



As listed in Table El the 35 C1/ 37 C1 abundance rati o is ra ther varied from around unity 



to around 3 and over. 



Cernicharo. Guelin. &: Kahane 



((2000) reported an X[ 35 C1/ 37 C1] ratio 



of 3.1±0.6 from line intensity ratios of NaCl, KC1, and A1C1 in IRC +10216. Recent H erschel 



HIFI observations show an X[ 35 C1/ 37 C11 ra tio of 2.7 in NGC6334I flLis et al. 112010( 1 and 



2.1 ±0.5 in W3A (ICernicharo et al. 



2010al ). These compare to the terrestrial value of ~3.1. 



It is interesting to note that the 35 C1/ 37 C1 abundance ratio comes in a good agreement with 
the integ ration line intensity ratio of the HC1 and H 37 C1 line, as predicted by the modeling 



efforts of 



Cernicharo et al. 



( )2010af ) for their Herschel HIFI observations of W3A. 



One may argue that the low isotopic ratio observed in W3IRS4 might be the result of 
observing an unresolved, very optically thick source. DR21(OH), which is at a comparable 
distance (~2.5 kpc vs. ~2.2 kpc for W3IRS4), with a comparable dust core size at 350/im 



(IDotson et al. 



20101 ) and HC1 opacity, argues against this scenario. 



Clearly, variations in HC1 isotopic ratio are highly localized. 



Salez. Frerking. fc Langer 



(119961 ) suggested one likely origin for such variations through explosive nucleosynthesis, 
whereby enrichment of one isotopologue over the other, as synthesized in massive stars, is 
mixed into the surrounding molecular cloud through explosive processes such supernova 



explosions. 



Models for nucleosynthesis in Type la supernovae (SIN 



Kobayashi et al. 



e la; 



Travaglio et al 

and references therein) and Type II supernovae (SNe II; (INomoto et al. 1 11997 



(120041) 



2006 



20061 )) indeed show a general deficiency of 37 CI in the processed material 
in most cases. In SNe la, the 35 C1/ 37 C1 abundance ratios varies from 3.5 to 5.4 among 
models of different ignition cond itions in the central white dwarf before the thermonuclear 



runaway (jTravaglio et al. 



20041 ). In SNe II, the elemental yields of the nucleosynthesis 



depend on the mass and metallicity of the progenitor, as well as the explosion energy. 
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Kobayashi et al. I (l2006t ) shows the 35 C1/ 37 C1 abundance ratios produced in SNe with high 
explosion energy (>10 E 51 with E 51 ~ 10 51 erg, so called hypernovae) are typically in the 
range of 1.1 — 2.8 for varying progenitor mass and metallicity, while the ratios varies in a 
larger range (to as high as 18) in normal SNe (with explosion enery ~Esi). It is inetresting 
to note that in normal SNe II with relatively large progenitor mass (>25M ) and metallicity 

abundance ratios are consistently below unity to as low as ~0.5 



(Z~0.02), the 35 C 


/ 37 C1 


(Kobavashi et al. 


2006) 



these ratios could support the observed range of the CI/ CI abundance ratio. 

W3IRS4, and possibly W3IRS5 present interesting exceptions of relative 37 C1 deficiency 
predicted in most nucleosynthesis models SNe. This likely confines the progenitors of the 
SNe to be stars of relatively high mass (>25M ) and high metallicity (Z~0.02). 



Kawabata et al. 



feoid ) and 



Perets et al. 



( 120 101 ) recently reported the discovery of a 



new category of subluminous supernovae. Although the origin for such faint supernovae is 
still controversial, these authors reported the unusual composition in the supernovae ejecta 
with rich helium and calcium. It would be interesting to examine the nucleosynthesis yields 
in this class of supernovae in regard of the 35 C1/ 37 C1 abundance ratios. 



5. CONCLUSION 

We conducted a comprehensive survey of HC1 J = 1 — line in the Galaxy. Of 
the 27 sources/positions observed, fourteen show emission, nine show absorption, two 
showed marginal detection in emission, and two are non-detections. Fourteen of the 
sources/positions were also observed in the H 37 C1 J = 1 — transition. 

RADEX simulations show that HC1 emission is mostly associated with warm dense gas, 
of order 10 6 cm -3 . HC1 abundance is fairly uniform, in the range 10~ 10 to a few xlO -9 , in 
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general agreement with pr e vious studies of OMC-1 by lSchilke. Phillips. &: Wangl (119951 ) and 
Salez. Frerking. fc Langerl ( 119961 ). However, 35 C1/ 37 C1 abundance ratios are found to vary 
considerably from around unity in W3IRS4 and W3IRS5 to ~5 in DR21(OH). They are in 
good agreement to the int egrated intensity ratios of the HC1 and H 37 C1 line as suggested by 



Cernicharo et al. 



(EOlOah 



The variations of the 35 C1/ 37 C1 abundance ratios are found to be highly localized. They 
could be supported by the varying yields of nucleosynthesis of supervoae from different 
progenitors. The low 35 C1/ 37 C1 abundance ratios seen in W3IRS4 and W3IRS5 likely 
confine the progenitors of the supernovae to stars of relatively high mass (>25M ) and high 
metallicity (Z~0.02). 



Caltech Submillimeter Observatory (CSO) is supported through NSF grant AST- 
0540882. 
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Fig. 1.— HCl and H 37 C1 J = 1 - spectra toward W3IRS4, W3IRS5, W3(OH), and two 
positions in OMC-1 offset from IRc2 in declination by -60"and -100", respectively. Velocity 
scale, Vlsr — Vo, refers to the main hyperfine component. 
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Fig. 2.— HC1 and H 37 C1 J = 1 - spectra toward MonR2, IRC+10216, IRAS16293, 
Sgr B2(M), and Sgr B2(N). Velocity scale, Vlsr — Vo, refers to the main hyperfme component. 
For sources with detectable continuum, baseline is not removed to show the continuum level. 
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Fig. 3.— HC1 and H 37 C1 J = 1 - spectra toward G5.89-0.39, W49N, W51FIR2, and 
DR21(OH). Velocity scale, Vlsr ~ Vo, refers to the main hyperfme component. For sources 
with detectable continuum, baseline is not removed to show the continuum level. 
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Fig. 4. — HC1 spectra toward rest of the sources in the sample. 
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Fig. 5. — A RADEX simulation run for OMC-1 (0,-60) with a fixed gas temperature of 45 K. 
Contours are the relative x 2 between the observed and modeled line temperatures, starting 
from 0.5 with a step of -0.1; blue and red curves mark the hyperfme component ratios under 
optically thin limit (3 between F=5/2-3/2 and F= 1/2-3/2 transitions in blue, 1.5 between 
F=5/2-3/2 and F=3/2-3/2 transitions in red). The solid smooth line and two dotted lines 
mark the observed optical depth of the HC1 line and its ±lcr error bars. 
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Table I. Source list of the survey 



Source 


RA(J2000) 


Dcc(J2000) 


V LSR 


HC1 


H 37 C1 


W31K,S4 


no/ioc ttl o ns n r\ 

U2 25 oU.9U 


^oOac/ on// on 

o2 Ud 2U. 29 


-42 


Y 


Y 


W31RS5 


02 25 40.46 


{ • <\ Of" F O 0^~7 

62 05 52.27 


-39 


Y 


Y 


Wj(OH) 


O O i\T n O OO 

02 27 03.82 


/~* "T P" O O A f~'f\ 

61 52 24.60 


-44 


Y 


Y 


OMC-1 


n p" o p" "io 

05 35 14.18 


or" oo on ft p" 

-05 22 30.55 


8 


Y 


Y 


OMC-1 (0, -60) 


n P" OP" "1 "IO 

05 35 14.18 


or" oo on p" p" 

-05 23 30.55 


8 


Y 


Y 


UMU-1 (U, -1UUJ 


n tr otr i/i io 

05 35 14.18 


ntr o/i io rrcr 

-U5 24 12.55 


8 


Y 


Y 


Orion Bar 


n tr opt n o nn 

05 35 23.99 


nrr o /i prn nT 

-05 24 59.97 


9 


Y 


N 


TV J" ^ , ^ T~> O 

MonHz 


r\n r\*7 An in 

06 07 46.19 


oo no on 

-06 23 08.30 


10.5 


Y 


Y 


Al 1 CjL9d1 


r\c O A OT A 

06 34 37.64 


r\ a 10 /iooo 

U4 12 4o.o2 


15 


Y 


N 


1RC+10216 


O o A r-T p- rr *") O 

09 47 57.30 


"IO 1 / 1 /lOOP" 

13 16 43.05 


-26 


Y 


Y 


CIT6 


10 16 02.34 


on • t a -1 o on 

30 34 18.30 


-2 


Y 


N 


IRAS16293 


16 32 22.91 


-24 28 35.52 


4 


Y 


Y 


Sgr B2(M) 


17 47 20.16 


-28 23 04.97 


62 


Y 


Y 


Sgr B2(N) 


17 47 19.90 


-28 22 17.78 


65 


Y 


Y 


Sgr A* 


17 45 40.04 


-29 00 28.10 


60 


Y 


N 


Gl.6-0.025 


17 49 20.61 


-27 34 05.46 


55 


Y 


N 


G5.89-0.39 


18 00 30.43 


-24 04 01.48 


10 


Y 


Y 


G10. 62-0.38 


18 10 28.67 


-19 55 50.08 


0.4 


Y 


N 


W43-MM2 


18 47 36.70 


-02 00 51.90 


90 


Y 


N 


W43-MM1 


18 47 47.00 


-01 54 28.00 


90 


Y 


N 


G34.26+0.15 


18 53 18.58 


01 14 58.41 


54.6 


Y 


N 


W49N 


19 10 13.62 


09 06 17.3 


0.0 


Y 


Y 


G45.12+0.13 


19 13 27.84 


10 53 36.6 


58.0 


Y 


N 


W51FIR2 


19 23 44.17 


14 30 33.4 


57.0 


Y 


Y 


DR21(OH) 


20 39 00.70 


42 22 46.7 


-4.7 


Y 


Y 



Table 2. Hyperfine structure fits for HC1 and H 37 C1 line emission. 



HCl H 37 C1 



Source 


VLSR 


w 


Si 


S 2 


s 3 


(AT-r) 


T 


Si 


S 2 


S 3 


(AT ■ r) 


T 


(AT)* 


r(HCl)* 


t(H 37 C1)* 


W3IRS4 


-44.0 


3.2 


4.0 


7.7 


5.9 


7.9(1.2) 


3.7(1.1) 


4.4 


8.0 


6.7 


6.7(0.6) 


1.8(0.6) 


2.7(0.8) 


2.5(2.1) 


2.8(2.5) 


W3IRS5 


-37.9 


4.8 


15.7 


24.9 


24.3 


24.6(1.9) 


5.8(0.7) 


12.8 


22.7 


17.8 


16.0(0.9) 


3.1(0.3) 


4.6(0.4) 


5.0(2.3) 


3.3(1.0) 


W3(OH) 


-49.2 


4.9 


7.9 


12.0 


12.2 


11.3(1.8) 


5.3(1.2) 


1.5 


9.2 


7.0 


3.9(0.5) 


0.7(0.7) 


2.2(0.6) 


4.9(6.9) 


1.9(2.3) 


OMC1 (0, -60) 


8.4 


3.1 


9.6 


17.6 


13.2 


15.1(1.4) 


2.0(0.5) 


4.9 


7.6 


4.4 


5.4(0.7) 


1.0(0.7) 


7.5(0.7) 


2.2(0.6) 


0.7(0.2) 


OMC1 (0, -100) 


7.9 


3.2 


10.7 


19.3 


17.1 


20.3(2.7) 


3.5(0.9) 


4.2 


8.5 


7.0 


6.0(0.6) 


0.7(0.6) 


6.1(0.7) 


3.6(1.6) 


1.0(0.3) 


ORIONBAR a 


10.6 


1.9 


1.2 


2.6 


1.8 


3.6(0.4) 


2.3(0.6) 












1.6(0.4) 


2.3(0.6) 




MonR2 


10.7 


1.7 


2.3 


4.8 


4.0 


9.0(0.6) 


3.3(0.4) 


1.1 


2.0 


2.0 


4.1(0.8) 


3.6(1.4) 


2.5(0.5) 


4.4(4.3) 


1.1(0.5) 


AFGL961 a 


13.4 


6.5 


1.3 


3.4 


2.3 


1.1(0.1) 


0.4(0.3) 












2.8(2.5) 


0.4(0.3) 




IRC+10216 


-17.8 


6.4 


5.2 


2.8 


5.5 


9.6(0.1) 


17.8(0.1) 


0.4 


0.8 


1.4 


0.8(0.3) 


5.0(3.0) 


0.54(0.01) 


7.3(6.8) 


0.8(0.5) 




-27.4 


8.5 


6.5 


5.9 


9.1 


6.0(0.1) 


6.5(0.1) 


2.0 


1.8 


2.0 


3.2(3.0) 


11.6(1.5) 


0.93(0.01) 


3.2(0.6) 


0.8(0.2) 


IRAS16293 


3.9 


4.2 


3.0 


8.2 


6.2 


4.5(0.3) 


1.1(0.4) 


1.4 


3.2 


2.3 


2.0(0.5) 


0.6(0.7) 


3.7(1.1) 


1.1(0.3) 


0.4(0.1) 


G5.89-0.39 b - d 


9.3 


4.9 


2.7 


8.0 


6.0 


3.5(0.7) 


0.7(0.5) 


1.0 


5.2 


3.9 


2.1(0.2) 


0.5(0.4) 


4.6(0.4) 


0.8(0.2) 


0.5(0.2) 




17.8 


13.5 


10.3 


17.0 


10.3 


2.8(0.2) 


0.3(0.2) 


3.2 


5.3 


2.4 


2.0(0.5) 


4.1(1.2) 


9.4(4.6) 


0.3(0.3) 


0.08(0.07) 




38.3 


5.4 


1.7 


5.3 


1.7 


1.6(0.2) 


0.1(0.1) 


0.6 


2.1 


0.4 


0.6(0.1) 


0.1(1.1) 


15(8) 


0.1(0.1) 


0.03(0.03) 


G5.89-0.39 c - d 


9.3 


4.9 


8.0 


24.0 


17.8 


10.8(2.1) 


0.7(0.5) 


2.8 


15.9 


11.7 


6.4(0.5) 


0.5(0.4) 


13.4(1.3) 


0.8(0.3) 


0.5(0.3) 




17.9 


13.3 


30.6 


50.6 


31.5 


8.7(0.5) 


0.4(0.2) 


9.7 


15.9 


7.2 


6.3(1.4) 


4.3(1.2) 


24(11) 


0.4(0.3) 


0.10(0.08) 




38.3 


5.4 


5.2 


15.8 


5.0 


4.8(0.3) 


0.1(0.1) 


2.0 


6.4 


1.1 


1.9(0.2) 


0.1(0.8) 


48(24) 


0.1(0.1) 


0.03(0.03) 


G45.12+0.13 a 


58.5 


3.8 


1.8 


2.8 


2.5 


3.1(0.7) 


5.2(1.8) 












0.6(0.2) 


5.2(1.8) 




DR21(OH) e 


-3.6 


4.7 


3.4 


6.9 


5.6 


4.4(0.4) 


2.6(0.6) 


0.6 


1.9 


0.9 


0.9(0.1) 


0.1(1.0) 


1.7(0.2) 


2.7(1.7) 


0.4(0.3) 



Note. — Numbers in parentheses are quoted error bars; Si, S2, and S3 are integrated intensity for the F = 1/2 — 3/2, 5/2 — 3/2, and 3/2 — 3/2 component of the 
J = 1 — transition, respectively. 

a (AT)* and r(HCl)* are from HFS fit of the HC1 line profile. 

b assuming no appreciable beam dilution. 

Accounting for a likely x3 beam dilution. 

d (AT)* for the 18.0 km s~ 1 and 38.3 km s _1 components are from the HC1 measurements alone as the low S/N in the H 37 C1 counterparts prevent reliable HFS fits. 
e (AT)* is from the HC1 measurement alone as the low S/N in the H 37 C1 counterparts prevent reliable HFS fit. 
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Table 3. RADEX fits for HC1 and H 37 C1 line emission. 



Source 




N(H 2 ) 


n(H 2 ) 


N(HC1) 


N(H 37 C1) 


X[HC1] 


X[ 35 C1/ 37 C1] 


K 


VV OlIXOl 


JO 


oy^o j 


(i o)(p) 






<H-±u; 


n s+ - 3 

U.8_ 2 


n q 


W3IRS5 


100 a 


5(23) a 


(8-13)(5) 


(1.5-2.2)(14) 


(1.3-1.8)(14) 


4(-10) 


i n+0.5 
L - Z -0A 


1.2 


W30H 


150 a 


4(23) a 


(2-3)(5) 


(1.5-2.4)(14) 


(7.0-9. 9)(13) 


5(-10) 


9 q+1.2 
z -°-0.8 


1.8 


OMC-1(0,-60) 


45 b 


7.4(22) c 


(l-2)(7) 


(5.1-6.3)(13) 


(1.8-2. 8)(13) 


8(-10) 


9 c+0.9 


2.4 


OMC-1(0,-100) 


30 b 


7.4(22) c 


(8-ll)(6) 


(1.3-1.5)(14) 


(5.5-7.8)(13) 


2(-9) 


9 i +0.5 


2.4 


Orion Bar 


75 d 


2.3(22) c 


(l-3)(6) 


(2.5-4.3)(13) 




l(-9) 






MonR2 


70 e 


3.9(22) e '^ 


(l-3)(6) 


(5.7-9.3)(13) 


(2.8-3. 7)(13) 


2(-9) 


9 q+1.0 
z -°-0.7 


2.2 


AFGL961 


20 9 


1.1(22)' 1 


(7-10)(6) 


(1.9-2.5)(13) 




2(-9) 






IRC+10216 


50 ! 


7(20)-* 


(5-30) (4) 


(5.7-8.8)(14) 


(2.6-3. 0)(14) 


l(-6) 


2 6 +as 


1.2 


IRAS16293 


so^ 


2(23) fc 


(l-4)(6) 


(3.0-4.7)(13) 


(1.3-2. 4)(13) 


2(-10) 


2 1+ 1 - 6 


2.5 


G5.89-0.39 1 


70' 


3.1(22)' 


(l-5)(6) 


(5.4-7.5)(13) 


(2.8-3. 3)(13) 


2(-9) 


9 i +0.6 


2.6 


G5.89-0.39 2 


70' 


3.1(22)' 


(3-13)(6) 


(7.4-10.1)(13) 


(3.2-3. 9)(13) 


3(-9) 


2 5+ - 7 


2.6 


G45.12+0.13 


25 m 


1.2(23) m 


(9-12)(5) 


(1.5-2.0)(14) 




l(-9) 






DR21(OH) 


40" 


3(23)" 


(l-3)(6) 


(7.1-11.3) (13) 


(1.2-2.5)(13) 


3(-10) 




4.7 



Note. — The numbers in the rightmost parentheses are exponents to the power of 10; R is the integrated intensity 
ratio of the HC1 and H 37 C1 J = 1 — transition. 

1 assuming no appreciable beam dilution. 

2 accounting for a likely x3 beam dilution. 
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